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Abstract 

We consider the following two deterministic inventory optimization problems over a finite planning horizon 
T with non-stationary demands. 

• Submodular Joint Replenishment Problem. This involves multiple item types and a single retailer who 
faces demands. In each time step, any subset of item-types can be ordered incurring a joint ordering 
cost which is submodular. Moreover, items can be held in inventory while incurring a holding cost. The 
objective is find a sequence of orders that satisfies all demands and minimizes the total ordering and 
holding costs. 

• Inventory Routing Problem. This involves a single depot that stocks items, and multiple retailer locations 
facing demands. In each time step, any subset of locations can be visited using a vehicle originating from 
the depot. There is also cost incurred for holding items at any retailer. The objective here is to satisfy all 
demands while minimizing the sum of routing and holding costs. 

We present a unified approach that yields O ^ lo^^iog t ) "f^'^tor approximation algorithms for both problems 
when the holding costs are polynomial functions. A special case is the classic linear holding cost model, wherein 
this is the first sub-logarithmic approximation ratio for either problem. 

Key words: Inventory Management, Approximation Algorithms, Submodular Function, Joint Replenishment 
Problem, Inventory Routing Problem 


1 Introduction 


Deterministic inventory theory provides streamlined optimization models that attempt to capture tradeoffs in man¬ 
aging the flow of goods through a supply chain. We consider two classical models in deterministic inventory 
theory: the Joint Replenishment Problem (JRP) and the Inventory Routing Problem (IRP). These inventory models 
have been studied extensively in the literature (see, e.g., |jl|], [|2T[]) and recently there has been significant progress 
on many variants of these models (see, e.g., [^, [22|, @], 0], [Oi). In this paper, we present a unified approach 


that yields approximation algorithms for both models with generalized cost structure - the JRP with submodular 
setup cost and the IRP with arbitrary embedding metric. 


The JRP with deterministic and non-stationary demand is a fundamental yet notoriously difficult problem in in¬ 
ventory management. In these models, there are multiple item types, and we need to coordinate a sequence of 
(joint) orders to satisfy the demands for different item types before their respective due dates. Ordering inventory 
in a time period results in setup costs (or fixed ordering costs), and holding inventory before it is due results in 
holding costs. The objective is to find a feasible ordering policy to satisfy every demand point on time over a 
finite planning horizon so as to minimize the sum of setup and holding costs. The JRP is a natural extension of the 
classical economic lot-sizing model that considers the optimal trade-off between setup costs and holding costs for 
a single item type (see [|^). With multiple item types, the JRP adds the possibility of saving costs via coordinated 
replenishment, a common phenomenon in supply chain management. 
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Most of the literature on deterministic JRP is on the additive joint setup cost structure. Under this structure, 
there is a one-time setup cost if any item type is ordered, and there is an individual item setup cost for each item 
type ordered; the joint setup cost for this particular order is simply the sum of the one-time setup cost and these 
individual item setup costs. The additive joint setup cost structure loses significant amount of modeling power 
and flexiblity (see p^, [p^, p^, p^]). In this paper, we adopt the joint setup cost structure introduced recently 
in [W] that satisfies fwo nafural properfies known as monotonicity and submodularity. The monofonify properfy 
means fhaf fhe joinf sefup cosf increases wifh fhe sef of ifem fypes ordered. The submodularify properfy capfures 
economies of scale in ordering more item fypes, i.e., fhe marginal cosf of adding any specific ifem fype fo a given 
order decreases in fhe sef of ifem fypes included. 


The IRP is also a classical problem in inventory managemenf fhaf capfures fhe frade-off befween fhe holding cosfs 
for invenfory and fhe roufing cosfs for replenishing fhe inventory af various locafions in a supply chain (see, e.g., 
0], lH], [0, @]). The problem involves multiple ifem fypes fhaf are stocked in a single depof, fhaf musf be 
shipped fo meef fhe demand for fhese ifem fypes arising af multiple refailers specified over fhe course of a planning 
horizon. Similar fo fhe JRP, fhe cosfs of holding a unif of invenfory af each refailer are specified fo compute fhe 
invenfory holding cosfs. Differenf fhan fhe JRP, we consider fransporfafion (or vehicle routing) cosfs in some mefric 
defined by fhe depof and refailers in fhe IRP, insfead of joinf sefup cosfs considered in fhe JRP. 


1.1 Main Results and Contributions 


We presenf a unified approach fhaf yields O ^ log log r ) -approximation algorifhms for bofh fhe JRP wifh submod- 
ular sefup cosfs and fhe IRP wifh any embedding mefric, when fhe holding cosfs are polynomial functions (which 
subsumes conventional linear cosfs as a special case). This is fhe firsf sub-logarifhmic approximafion rafio for 
eifher problem under fhese cosf sfrucfures. 

We remark fhaf if fhe sefup cosf function in submodular-JRP is fime-dependenf fhen fhe problem (even wifh zero 
holding cosfs) becomes as hard fo approximate as sef cover [01. The same observation is frue if fhe mefric in IRP 
is fime-dependenf. So our sub-logarifhmic rafio approximafion algorifhm relies crucially on fhe uniformify of fhese 
cosfs over time. 


For fhe submodular JRP, Cheung ef al. [10] obfained consfanf-facfor approximafion algorifhms under several spe¬ 
cial submodular funcfions (i.e., free, laminar and cardinality). In confrasf, we consider general submodular func¬ 
tions wifh special (polynomial) holding cosfs. 


For fhe IRP, Fukunaga el al. [ jT^ considered a restricted class of “periodic policies” and obtained a constant-factor 
approximation algorithm. Whereas our result is for arbitrary policies and polynomial holding costs. 


A straightforward modification of our algorithm for polynomial holding costs also yields 0(log T)-approximation 
algorithms for submodular JRP and IRP with arbitrary (monotone) holding costs. The submodular JRP result 
improves upon the approximation ratio of 0{log{NT)) by [p^. The IRP result is incomparable to the O(logri) 
approximation ratio mentioned in [|l7|], where n is the number of retailers. 


1.2 Our Approach 

At a high-level, the algorithm for submodular JRP has the following steps. (The algorithm for IRP is very similar 
- we in fact present an algorithm for a unified problem formulation.) Firsf, we solve a nafural lime-indexed LP 
relaxalion fhaf was also used in [0]. Then we consfrucf a “shadow interval” for each demand poinf fhaf corresponds 
fo fractionally ordering half unif of fhe ifem. We also sfrefch each shadow interval appropriately (depending on 
fhe degree of fhe holding cosf funcfion) so as fo obfain an oplimal Irade-off befween holding and sefup cosfs: Ihis 
is whal resulfs in fhe O ^ lo^g fog t ) approximation rafio. Nexf, we parlifion fhese sfrelched intervals info multiple 
groups based on well-separated widfhs. Finally we place a separafe sequence of orders for each group, and argue 
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using submodularity of the setup cost function that the total setup cost of each group is bounded by the LP setup 
cost. This step relies on the fractional subadditivity property of submodular functions. 


It turns out that we do not require the full strength of submodular functions: the algorithm and analysis work even 
for functions satisfying an approximate notion of fractional subadditivity (see Definition ^ ) as long as the natural 
LP relaxation can be solved approximately. This allows us to also obtain an approximation algorithm for IRP 
since the TSP cost function satisfies 1.5-approximafe fracfional subaddifivify and fhere is a 2 + o(l) approximafion 
algorifhm for ifs LP relaxafion (see Secfion ^ for defails). 


We believe fhaf some of our techniques may be useful in obtaining a constant factor approximation algorithm for 
both problems in their full generality. 


1.3 Literature review 


As mentioned earlier, most of the existing literature on deterministic JRP with non-stationary demand uses the 
additive joint setup cost structure. Arkin et al. [ffl] showed that the additive JRP is NP-hard. Nonner and Sviridenko 


[ 24 ] further showed that the additive JRP is in fact APX-hard with nonlinear holding cost structure. There have 


been several approximation algorithms for the additive JRP (see | |10| ] and the references therein). The state-of-the- 
art approximation algorithms for the additive JRP are due to [^], [^] and |^], with approximation ratios 2, 1.80 
and 1.791, respectively. 


Due to the limited modeling power of the additive JRP, Cheung et al. [ ]10| ] first studied the submodular JRP in 
which the setup costs are submodular. They gave an 0(log(A^T))-approximation algorithm for the general sub¬ 
modular JRP (where N is number of items and T is number of periods). They also analyzed three special cases 
of submodular functions which are laminar, tree and cardinality cases. They showed that the laminar case can be 
solved optimally in polynomial time using dynamic programming, and obtained a 3-approximation for the tree 
case and a 5-approximation for the cardinality case. Our work contributes to the literature by giving approximation 
algorithms for the general submodular JRP with special holding cost structures. 

The IRP has also been studied extensively in the literature (see [i> [[T^], [[T^, [ |T7| ] for an overview of this problem). 
The problem can be cast a mathematical program (see, e.g., ||^) and solution approaches typically involve heuris¬ 
tics that trade-off between holding and transportation costs (see 10, |], 0, l[n]], [0, 0). Closer to our work, 
Fukunaga et al. [ [T^ gave constant factor approximation algorithms for the IRP restricting to periodic schedules. 
In contrast, our results do not require the schedule to be periodic but require polynomial holding costs. 


1.4 Structure of this paper and some notations 

We organize the remainder of the paper as follows. In Section we present a unified formulation for fhe submodu¬ 
lar JRP and fhe IRP wifh arbifrary embedding mefric, and sfafe our main resulf. In Secfion^ we propose a unified 
approximafion algorifhm for bofh problems. In Section ^ we discuss how fo solve fhe LP relaxafion efficienfly. 
We conclude our paper in Section |5[ 

Thr'oughouf fhe paper, we use fhe nofafion [xj and [x] frequenlly, where [xj is defined as fhe largest integer 
value which is smaller than or equal to x; and [x] is defined as fhe smallesf infeger value which is greater fhan or 
equal fo x. Additionally, for any real numbers x and y, we denote x"*" = max{x, 0}, x V y = max{x, y}, and 
X Ay = min{x, y}. The nofafion := reads “is defined as”. 


2 A Unified Formulation for the JRP and the IRP 

In fhis secfion, we formally describe a unified problem sfafemenf fhaf includes fwo classical deterministic inven¬ 
tory problems as special cases, i.e., fhe joinf replenishmenf problem (JRP) wifh submodular sefup cosfs and fhe 
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inventory routing problem (IRP) with arbitrary embedding metric. We also present a unified framework for this 
problem, and state our main result. 


2.1 Problem Statement 


There are N elements (e.g., item types in the JRP or retailers in the IRP) that are needed to serve external demands 
over a finite planning horizon of T periods; these elements are denoted by the ground set Af = N}, and 

the time periods are denoted by the set T = {1,... , T}. For each time period t £ T and each element i £ Af , 
there is a known demand da > 0 units of that element. We use V to denote the set of all strictly positive demand 
points (f, t) with da > 0. To satisfy these demands, an order may be placed in each time period. Each demand 
point (f, t) £ V has to be served by an order containing element i before or at time period t, i.e., no backlogging 
or lost-sales are allowed. 

The inventory system incurs two types of cost - the joint ordering cost and the holding cost. 


• The joint ordering cost is a function of the elements that place strictly positive orders in any given period. 
More specifically, for any time period t and a subsef of elements S C J\f, the joint ordering cost of ordering 
demand for elements in S in period f is a function of S, which is denoted by f{S). 

• Because the setup cost of ordering an element is independent of the number of units ordered, there is an 
incentive to place large orders to meet the demand not just for the current time period, but for subsequent 
time periods as well. This is balanced by a cost incurred by holding inventory over time periods. We use /i*^ 
to denote the holding cost incurred by ordering one unit of inventory in period s, and using it to meet the 
demand for element i in period t. We assume that is non-negative and, for each demand point (f, t), is a 
nonincreasing function of s, i.e., holding inventory longer is always more costly. Thus, if the demand point 
(z, t) is served by an order at time period s, then the system incurs a holding cost of Hl^ := dahl^. 


The goal is to coordinate a sequence of (joint) orders to satisfy all the demand points on time so as to minimize the 
sum of joint ordering and holding costs over the T periods. 

The above unified problem stafemenf encompasses fwo classical deterministic inventory problems described below. 


The submodular JRP. The JRP involves multiple item types and a single retailer who faces demands. In each 
time step, any subset of item-types can be ordered incurring a joint ordering cost which is submodular. The 
objective is find a sequence of orders thaf satisfies all demands and minimizes the total ordering and holding costs. 
The elements in the above problem statement are the item types in the JRP; and the joint ordering cost /(•) is 
commonly referred to as the setup cost (or equivalently, the fixed ordering cost) in the JRP. 


The submodular JRP considers a special class of /(•) called submodular functions (see, e.g., [10]). More precisely, 
we assume that the function /(•) is non-negative, monotone non-decreasing, and also submodular. The non¬ 
negativity and monotonicity assert that for every 5i C 52 C J\f, we have 0 < f{Si) < f{S 2 )- Submodularity 
requires that for every set Si, S 2 ^ M, we have 


f{Si) + f{S 2 ) > f{Si u 52) + /( 5 i n 52). 


There is an equivalent definition that conveys the economies of scale more clearly. That is, for every set Si C 
S 2 and any item type i £ M,we have f{S 2 U {z}) — f{S 2 ) < f{Si U {z}) — f{Si), i.e., the additional cost 
of adding an item type to the joint order is decreasing as more item types have been included in that order. 

The IRP with arbitrary embedding metric. The IRP involves a single depot r that stocks items, and a set of 
retailer locations (denoted by the ground set M) facing demands. In each time step, any subset of locations can be 
visited using a vehicle originating from the depot. The objective here is to satisfy all demands while minimizing 
the sum of routing and holding costs. The elements in the above unified problem statement are the retailers in the 
IRP; and the joint ordering cost /(•) is the shipping or routing cost. 
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The IRP is specified by a complete graph on vertices V with a metric distance function M+ that satisfies 

symmetry (i.e. w{ba) = w{ab) for any a,b ^ V) and triangle inequality (i.e. w{ab) + w{bc) > w{ac) for any 
a,b,c £ V). The vertex set 1/ = AA U r, containing the depot and the set of retailers. The shipping or routing cost 
f{S) can be defined as the travelling salesman (TSP) cost of visiting the retailers in 5 C 1/. Formally, 

f{S) := minimum length of tour that visits each vertex in 5 U {r}, MS^M. (1) 


2.2 IP Formulation and its LP Relaxation 


The unified problem described above can be written as an integer programming problem as follows (see also [HD- 
First we define two types of binary variables yf and such that 




1, if the subset of elements 5 C is ordered in period s, 

0, otherwise. 

1, if the demand point (t, t) is satisfied using an order from period s, 
0, otherwise. 


Then the integer programming (IP) formulation is given by 


(IP) 


T t 

min 


SCA/" s=l 


Y^'st = D 

S = 1 

V(i,f) G V 

x\t < Y y 

f, V(z,f) GP,Vs = l,...,f 

S:ieScN 


A Vs e 

V(z,t) gT>,Vs = 


( 2 ) 


The first constraint in (^ enforces that every demand point {i, t) must be served by an order before or in time 
period t. The second constraint in (^ ensures that the joint order S has to contain element i if any demand (z, t) 
is served at time period s. There is a natural linear programming (LP) relaxation of (IP) that relaxes the integer 
constraints on xl^ and yf to non-negativity constraints. 

To obtain approximation algorithms for the IP (^ using our framework, we only need to assume that the set function 
/(•) satisfies an approximate notion of fractional subaddivitity (which is much weaker than submodularity). 

Definition 2.1 (/3-approximate fractional subadditivity) The set function /(•) is /3-approxiniate fractional sub¬ 
additive, if for any collection {Si, Aj} of weighted subsets with 0 < A* < 1 and Yli\veSi v £ S, 

we have f{S) < /3 • ^ Xif{Si). Namely, if the sets Si form a fractional cover of S f M, then the cost of S is at 
most (3 times the sum of the costs f{Si) weighted by the corresponding coefficients. 


It is known that if a function is submodular, then it is also fractional subadditive (see [16]), i.e., the notion of 
submodularity is stronger. For the submodular JRP, the setup cost function f{S) is submodular and hence also 
fractional subadditive (or equivalently, 1-approximate fractional subadditive). 


For the IRP with arbitrary embedding metric, the vehicle routing cost f{S) although not submodular, can be shown 
to be 1.5-approximate fractional subadditive. This follows from the fact that the natural LP-relaxation for TSP has 
an integrality gap at most 1.5 (see [^] and [^]). 


Note that the LP relaxation of (^ has an exponential number of variables; we need to ensure that this LP relaxation 
can be (at least approximately) solved efficiently. 
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Definition 2.2 ( 7 -approximate LP solution) Wfe say that a feasible LP solution is y-approximate, if its objective 
value is at most 7 times the optimal LP objective value. 

Using the ellipsoid method one can compute efficiently: an exact LP solution for the submodular JRP and a 
(2 + o(l))-approximate LP solution for IRP We delegate the discussion on this to Section ^ for better readability 
of this paper. 


2.3 Our Main Results 


Assumption 2.3 (a-degree polynomial holding cost) For each element i £ Af and 1 < s < t £ T, the holding 
cost of holding an inventory unit of element ifrom period s to t is 

Kt = [t-sTK 


for some base per-unit holding cost h\> and some a > 1 . 


Note that when a = 1, this reduces to the conventional linear holding cost. We also have = d\{t — s)°^h\. 
Now we are in a position to formally state our main result (which will be proved in the following section). 

Theorem 2.4 Under Assumption |2. 4 there is an O -approximation algorithm for the integer pro¬ 

gram defined in provided that /(•) is ^-approximate fractional subadditive, and a y-approximate solution to 
the LP relaxation of can be found in polynomial time. 


Corollary ^ below is an immediate consequence of Theorem ^ , since 


1. /3 = 7 = 1 for the submodular JRP; 

1. (5 = 1.5 and 7 = 2 + o(l) for the IRP. 

Corollary 2.5 Under Assumption there is an O -approximation algorithm for both the submodular 

JRP and the IRP with arbitrary embedding metric. 


To the best of our knowledge, this is the first sub-logarithmic approximation ratio for either problem. 

We remark that it is immediate that our approach yields O (log T)-approximation algorithms for submodular JRP 
and IRP with arbitrary (monotone) holding costs (i.e., waiving Assumption 2.3). 


3 LP-Rounding Algorithm 


We present an LP-rounding algorithm for the integer program (g) under Assumption in Section |3.1| , and then 
carry out a worst-case performance analysis in Section 3.2. 


3.1 Algorithm Description 

We describe our procedure of rounding a 7 -approximate solution (y, x) of (LP). We set p := [(log . 
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Step 1 - Constructing extended shadow intervals. We first construct what-we-call extended shadow interx’als 
as follows. For each demand point {i, t), we take the 7 -approximate LP solution and find a time period such 
that 

t t 


x\t > 1/2 and 


E < 1/2, 

o' _L1 


i.e., finding the closest s to the left of t such that the sum of x variables for (z, t) contains a half point. Then , t] 
is called the shadow interval for this particular demand point (z, t). We also measure its length := t — s'^- 

Next for each demand point (z, t), we round the length up to the nearest power of p. If t) ~ ^ 

s*- y = t. Else we find the smallest integer m>l such that >t — s'y ^y and then stretch the original shadow 

interval from t to where 

(2,t) 

We call the interval [s*- ^yt] the extended shadow interval for the demand point (z, t), and also measure its length 
Figure || below gives a graphical representation of this step. 


first time sum of x'^t crosses 1/2 demand d't>0 

i 1 

I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I 

S*(i,t) ^ 


extending the length of shadow interval 
for (i,t) to the nearest power of p 


Figure 1: Illustration of an extended shadow interval for demand point (z, t). 


Step 2 - Partitioning demand points according to extended shadow intervals. Next we partition the demand 
points according to the length of their extended shadow intervals. For each demand point (z, t), its length I*- y falls 
into exactly one of the values below (recall by construction I*, y is either zero or an integer power of p). 

AT}, where A: = 1+ [log^T] . 

In this way, we have partitioned the demand points into k = O ^CT io^gfogy ^ number of groups as follows, 

£0 = {(At) = 0 } and = {(z, t) G P : ^ = p™ A T)} , V m G {1,... , A: - 1}. 

That is, the shadow intervals within each group Cm share the same length: 

r 0 if ZTT, = 0 

:= I A T if m G {1, 2,... , A: - 1} 

Step 3 - Placing orders. Based on the above partition of demand points, we describe our ordering procedure. 

Now fix an m G {0,1,..., A: — 1} and focus on the demand group Cm- Fet Tj = 1+j ■ Wm (< T) for j = 0,1,.... 
We place a tentative (joint) order in each period Tj (j = 0,1,...), i.e., once every Wm periods. In each period 
Tj < T (j = 0,1,2 ...), we identify the set of elements 

= {* : (At) G Cm and Tj G [s(i,t),t]} , 
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i.e., all the elements within Cm whose shadow intervals contain (or intersect with) time period Tj. We then place 
an actual joint order that serves the demand points associated with in period Tj. Figure ^ gives one specific 
example of how the algorithm places these orders. 

We repeat the above procedure for all groups m = — 1. In any given period, if there is more than one 

joint order (across different groups), we simply merge them into a single joint order. 

Focusing on the set of demand points Eiements: item types in JRP / 

retaiiers in the iRP 


shadow intervais: 


time periods: 


eiement 1 

element 2 element 4 

element 4 elements 

element 2 

I_I_I_I_I_I_I_I 

1 Ti T2 Ts T4 Ts Ts ... T 


elements ordered: {2,4} 


{1} {4,3,2} 


Figure 2: Placing actual orders for the demand points within Cm 


This concludes the description of our LP-rounding algorithm. 

3.2 Worst-case Analysis 

We shall prove that our LP-rounding algorithm gives an O -approximation for the unified problem. 

For brevify, we jusl call fhe y-approximafe solution (y, x) fo fhe LP relaxafion of (|^ the ^-approximate LP solution. 


Analysis of Holding Cost 

Lemma 3.1 Let ^yt\ be the shadow interval for some demand point (i, t). Then we have 

t 

{i,t) 


(3) 


S = 1 


where ’s are in the ^-approximate LP solution. 

Proof: By fhe consfrucfion of shadow infervals, we have Yl\=s' 




<1/2. Now since = 1 by the 


firsf consfrainf in (^, we musf have > 1/2. Hence we have 








S = 1 


I I 

I f 
(i,ty 


5 = 1 


V 


it<i I = 


2-E^ 


* A 


(4) 


s=i 




5 = 1 


where fhe second inequalify is due fo < Hgt for all s < s'^. by monofonicify of holding cosfs. 

Lemma 3.2 The total holding cost for the solution found by the LP-rounding algorithm is at most 

t 

0(0^)- E E<.<. 


where x\^ ’s are in the y-approximate LP solution. 












Proof: By the polynomial holding cost structure, for each demand point (i, t) in some C^, we have 


HI*. = diit-sy^,^r = diiwmrh\ 

HI, = diit-s[,^,^r>diiwm/prh\ 

(i,t) y ’ ^ 

where the inequality follows from the construction of extended shadow intervals. Hence it is clear that 





(i,i) 


t' 


By the LP-rounding algorithm, for each demand point we must have placed a (joint) order containing element 
i inside its extended shadow interval t)H]- Due to monotonicity of holding costs, the worst-case (that gives the 
highest possible holding cost) happens when our algorithm places the order at exactly time period s*- to satisfy 
the demand point (i, t). Hence, the total holding cost associated with the demand point {i, t) is upper bounded by 


H 


ii,t) ^ ® 

S = 1 


sHst^ 


where the second inequality follows from Lemma |3.1|. Now setting p 


[(log yields the result. ■ 


The intuition behind Lemma ^[2| is that when we stretch the shadow interval to the nearest integer power of p, 
the holding cost within the extended shadow interval does not grow too large due to the polynomial holding cost 
structure. In particular it grows by at most a factor of O(yTogT). On the other hand, stretching the shadow 
intervals in this manner gives us a tighter bound on the ordering cost (as shown below). 


Analysis of Ordering Cost 

To analyze the ordering cost component, we introduce the following bridging problem: 


(Covering-LP) min 


s.t. 


E E 

s<^j\r s=i 
t 

Y ^ y{i,t)ev 

zf > 0, 


(5) 


Vs = l,...,f,VS’CA(. 

The intuition behind introducing this bridging problem is as follows: if our algorithm places an order to satisfy 
the demand within its extended shadow interval then Lemma 3.2 implies that the holding cost can be bounded by 
0(-^log T) times the LP holding cost. Thus, the problem reduces to finding a “cover” for these intervals as defined 
in Problem (Q). In fhe remainder of fhe worsf-case analysis, we will focus on analyzing Ibis Covering-LP. 

Lemma 3.3 The optimal objective value of the Covering-LP is at most 

T 


2 - EE/(®te 

scat s=l 


s 


where yf ’s are in the -approximate LP solution. 
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Proof: We first check zf = (where yf is the 7 -approximate LP solution) is feasible to the Covering-LP 
defined in (^. It is obvious that zf = 2yf > 0 for all s = 1,... , t and for all S C Af. It suffices to verify the first 
set of constraints. Indeed, for each {i,t) G V, we have 

E E E E 29f> E 24i>i. (6) 

S:i€SCAr S-ieSCAf 

where the first inequality follows from the second constraint in (§), and the second inequality follows from the fact 

that Yll=si. , > Yll=s' . > 1/2 (by the definition of shadow intervals and their extensions). 

(i.t) (i,t) 

Hence, the optimal objective value of the Covering-LP 

f{s)zi <y.y. = 2 E E 

SQM s=l SQM s=l SQM s=l 

where the first inequality follows from that zf is optimal while zf is feasible. ■ 

Fix an m G {0,1,... , fe — 1} and we focus our attention on the demand group Cm (which have equal length 
of shadow intervals). We shall show that the total ordering cost associated with the set Cm by our LP-rounding 
algorithm can be upper bounded by 2/3 times the Covering-LP cost. Our proof strategy relies on the notion of 
approximate fractional subadditivity (see Definition ^Tl] ). 

Lemma 3.4 The total ordering cost associated with the set Cm by our LP-rounding algorithm is at most 

2/3- EE/(^)^^ 

SQM s=l 


where zf ’s are the optimal Covering-LP solutions. 

Proof: Recall that for each demand group Cm, the LP-rounding algorithm places a tentative (joint) order in each 
period Tj < T [j = 0, 1,...). Then in each time period Tj the algorithm identifies the elements Am within Cm 
whose shadow intervals contain (or intersect with) Tj and places an actual (joint) order Am that includes all of 
these “intersecting” elements. 

Now take any Tj < T: since the length of the shadow intervals in Cm is exactly Wm, all the shadow intervals 
associated with the order Am must lie within the interval (rj_i,Tj+i] (see Figure ^ as an example). Our LP- 
rounding algorithm places an actual (joint) order Am in period Tj and incurs an ordering cost f[Am)- We will 
show that the Covering-LP provides us with a fractional cover of Am which will be used to upper bound f{Am)- 

Indeed, for each demand point (i, t) associated with the order Am, we have 

"G'+i t 

E E E H d> t. E d>r (8) 

S:i€SCMs>rj.l s>Tj.i S:ieSCM s=s'. S:ieSCM 

- J J - _ 

where the first inequality holds because every shadow interval associated with the order Adm must lie within the 
interval (Tj_i, Tj+i] and the last inequality follows from the first constraint in the Covering-LP (^. 

Since /(•) is ^^-approximate fractional subadditive, then according to Definition |2.ij 

'^7 + 1 

E E (9) 

S:SQM s>Tj-i 
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It is then immediate that the ordering cost associated with the set Cm by our LP-rounding algorithm 


Tj + l 


j>0 


f(AL) < /S' E E E 


S:SCJ\f j>0 s>Tj-i 


z^fiS) < 2 / 3 - 

S:SCJ\rs=l 


(10) 


Lemma 3.5 The total ordering cost for the solution by our LP-rounding algorithm is at most 

T 


O al3 


^ogT 
log log T 




scj\rs=i 


where yf ’s are in the 'y-approximate LP solution. 


Proof: By Lemmas ^ and for each group Cm (m = 0,1,... , /c — 1), we conclude that the total ordering cost 
associated with the set Cm in our LP-rounding algorithm is at most 

T T 

2/3 • E E < 4/3 • E E 

SCAT S=1 SCAT S=1 


where Zg ’s are the optimal Covering-LP solution and y* ’s are the y-approximate LP solution. Then the result 

IorT \ 


follows from the fact that the number of groups k = O ( ) ■ 


Now we are ready to prove our main result Theorem 2.4. 


the solution by our LP-rounding algorithm is at most O 


Proof of Theorem Combining the results from Lemmas 3.2 and 3.5, the total holding and ordering costs for 

' lo^gioE ) times the 7 -approximate LP solution. ■ 


Remark: The O ^ iog°ioE ) approximation ratio is the best tradeoff achievable (in our approach) between the loss 
in holding and ordering costs, even under linear holding costs. Recall that for a given set W of widths for extended 
shadow intervals, the loss in ordering cost is just the number |VP| of distinct widths and the loss in holding cost 
depends on the aggregate stretch-factor incurred when the width of each shadow interval is increased to a value in 
W. Even if we allow for an arbitrary set W of widths (that may depend on the LP solution) and compute the worst 


ratio (using a “factor revealing linear program” as in [Il9|]) then we obtain O f ) as the approximation ratio. 


A Special Case with Perishable Goods 


We now consider a special holding cost which models perishable items with a fixed life-time c > 0. For each 
demand point (z, t), we can only start satisfying this order c periods before t, i.e., the ordering window is [t — c,t]. 
This setting is equivalent to the following holding cost structure. For each i £ Af and 1 < s < f G T, 


i _ 

st 



if t — c < s < f, 
if s < f — c. 


We also have Hl^ = dlhl^. 

In this setting, for each demand point (i, t), the extended shadow interval is simply [t — c, t] with length c. Hence 
our LP-rounding algorithm and its worst-analysis will apply with just a single group, and we obtain: 

Theorem 3.6 When items are perishable with a fixed life-time and the holding cost is negligible, the LP-rounding 
algorithm gives a 2-approximation for the submodular JRP, and a (6 -|- o{l))-approximation for the IRP 
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4 Solving the LP Relaxation 


As mentioned earlier in Section ^ the LP relaxation has an exponential number of variables and we first argue that 
there is an efficient way of solving this LP. We can readily write the dual of (LP) as 

(DLP) max ^ (11) 

s.t. h\<Hl, + bi,, V(f,f) GP,Vs = l,...,f 

T 

Vs = i,...,r,vsc iv} 

ies t=s 

5*t>o, y{i,t)€V,ys = i,...,t. 

The bl and 5*^ are the dual variable corresponding to the first and second constraints in the LP relaxation of (^). 
Note that the dual formulation (pT|) has an exponential number of constraints. 


Submodular JRP 

In the submodular JRP, the left hand side of the second constraint f{S) — Xligs Ylt=s Kt clearly submodular. 
Thus, there is an efficient separation oracle by using submodular function minimization to find violafed con- 
sfrainfs. This implies thaf fhe dual problem (and therefore the primal) can be efficiently solved using the ellipsoid 
method. This was also discussed in [|T^]. 


Approximately Solving the LP for IRP 


The TSP costs /(•) are not submodular, and in fact the above separation problem is NP-hard. However, there is an 
approximate separation oracle (see Lemma 4.3) which suffices fo compufe an approximafely optimal solution to 
(DLP). The main ingredient is an approximation algorithm for the following auxiliary problem: 


Definition 4.1 (Minimum ratio TSP) The input is a metric {V, w) with a designated depot r G IV and rewards 
a : L —)• M+ on vertices. The goal is to find a subset S 'TV that minimizes: 

where f{S) is the TSP cost as defined in (^anda{S) = '^ai. 

ies 

There is a (2 + o{l))-approximation algorithm for the minimum ratio TSP problem. 


m 

a(5)’ 

Theorem 4.2 (Garg iJUl) 


Proof: The algorithm for minimum ratio TSP uses the 2-approximation algorithm for the related A:-TSP problem 
(i.e. given a metric, depot r and target k, find a minimum lengfh four from fhe depof fhaf visifs af leasf k vertices). 
The algorifhm which is based on sfandard scaling argumenfs, is given below for complefeness: 


1. Guess (by enumerating over \ V\ choices) the maximum reward vertex u in an optimal solution. 

2. Remove all vertices with reward more than a„. 

3. For each n G IV set its new reward d^ to be the largest integer such that • ^ < a„. 

4. For each k = 1, • • • ,n^, run the /c-TSP algorithm with target k on the modified metric containing dv co¬ 
located vertices at each v £ V. 

5. Output the best ratio solution found (over all choices of u and k). 


12 




It is easy to see that this algorithm runs in polynomial time since each < in?. We now show that it has an 
approximation ratio of 7 = 2 + o(l). Let S* denote an optimal solution and u £ S* the maximum reward vertex. 
Consider the run of the above algorithm for this choice of u: note that none of the vertices from S* is removed. By 
the definition of new rewards, we have a„ — ^ < • ^ < for all v G S*. So ^ ~ ^ — 

(1 — l/n)a{S*), which implies (as a is integer valued) that d(5*) > k := [^(1 — For this choice of 

k, the fe-TSP algorithm is guaranteed to find a subsef S Q V wifh a{S) > k and f{S) < 2 • f{S*). The rafio of 
fhis solution is: 

/(5) ^ 2/(5*) 2/(5) ^ 2 fjS*) 

a{S) ~ a{S) ~ Qu d(5) “ 1 — 1/n a{S*) 

Hence fhe above algorifhm achieves a 2 + o(l) approximafion guaranfee for minimum rafio TSP. ■ 


Lef 7 = 2 + 0 ( 1 ) denofe fhe approximafion guaranfee for minimum rafio TSP from Theorem O. We now show 
fhaf fhis algorifhm leads fo an approximafe separation oracle. 


Lemma 4.3 There is a polynomial time algorithm, that given vectors b = {b\ : (i, t) G P} and b = {h\^ : (i, t) G 
P, 1 < s < f}, outputs one of the following: 


1. A constraint in (DTP) that is violated by (b, b). 

2 . Certificate that (7b, 7b) is feasible in (DTP). 


Proof: The firsf sef of consfrainfs in (DLP) and non-negafivify of are easy fo verify since fhey are polynomial 
in number. Below we assume fhaf fhese are safisfied. 


In order fo verify fhe second sef of consfrainfs in (DLP), we use Theorem p~2| . For each s G [T] define an insfance 
Xg of minimum rafio TSP on mefric (V, w), depof r and wifh rewards '■= Ylt=s ^st for v £ V. Lef Ag denofe 
fhe solufion found by fhe approximafion algorifhm of Theorem 1.2 on Xg. 


If any solufion Ag has rafio less fhan one fhen if provides a violated consfrainf for (b, b). This corresponds fo fhe 
firsl condition in fhe lemma. 


If every solufion Ag has rafio af mosf one, we will show fhaf fhe second condifion in fhe lemma holds. The non- 
negafivify consfrainfs are clearly safisfied by (7b, 7b). To check fhe firsf sef of consfrainfs in (DLP), nofe fhaf for 
any {i, t) £V and s £ [T], 

^ {b\ - b\t) < max { 0 , b\ - 6 * J < 


To check fhe second sef of consfrainfs, nofe fhaf for any s £ [T] we have by fhe approximafion guaranfee in 

1 


Theorem 4.2, 


mm — 


f{S) 


> 


ves Ya=s Kt 


This implies fhaf (^b, ^b) satisfies all fhese consfrainfs. 


Using fhe above separafion oracle for (DLP) wifhin fhe ellipsoid algorifhm, we obfain a 7 -approximafely opfimal 
solufion fo (DLP), see eg. ||20|]. Then solving (LP) resfricfed fo fhe (polynomially many) variables fhaf are dual fo 
fhe consfrainfs generafed in solving (DLP), we obfain a 7 -approximafely opfimal solufion fo (LP) as well. 


Running time. Using fhe linear programming algorifhms in ^ along wifh some preprocessing, fhe running 
fime of fhe above approach is dominated by 0{V^T^) plus 0{T>T‘^) calls[] to a subroutine for: 

• submodular function minimization in case of JRP 

• minimum-ratio TSP in case of IRP. 

’The O notation hides logarithmic factors. 
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5 Concluding Remarks 


We presented an O ^ lo'giogT ) -approximation algorithm for submodular-JRP and IRP when holding costs are poly¬ 
nomial functions. Moreover, this approach applies to any ordering cost for which the correponding LP relaxation 
can be solved approximately and the ordering cost satisfies an approximate notion of fractional subadditivity. Ob¬ 
taining a constant-factor approximation algorithm for submodular JRP and IRP on general metrics (even with linear 
holding costs) remain the main open questions. 
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